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Abstract 

We analyze the critical transitions (a) to endemic states in an SIS epidemiological 
model, and (b) to full synchronization in an ensemble of coupled chaotic maps, on 
networks where, at any given time, each node is connected to just one neighbour. In 
these "monogamous" populations, the lack of connectivity in the instantaneous in- 
teraction pattern -that would prevent both the propagation of an infection and the 
collective entrainment into synchronization- is compensated by occasional random 
reconnections which recombine interacting couples by exchanging their partners. 
The transitions to endemic states and to synchronization are recovered if the re- 
combination rate is sufficiently large, thus giving rise to a bifurcation as this rate 
varies. We study this new critical phenomenon both analytically and numerically. 
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1 Introduction 



The spontaneous emergence of different kinds of collective behaviour is the 
most paradigmatic phenomenon in natural and artificial systems formed by 
large ensembles of interacting dynamical elements. The interplay of individual 
dynamics and interactions entrains elements into coherent macroscopic evo- 
lution, which typically manifests itself in the form of spatial and/or temporal 
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structures. Pattern formation and synchronization are widespread examples 
[1,2]. 



Self-organization into collective evolution requires that the information on the 
state of any single element be able to spread all over the system, eventually 
reaching any other element in the ensemble. A crucial ingredient that controls 
such mutual influence of any pair of elements is the interaction pattern of 
the ensemble. In the case of binary interactions, this pattern is conveniently 
represented as a network, whose links join pairs of elements which interact 
with each other [3] . Coherent evolution of the whole ensemble is possible only 
if the interaction network is not disconnected. 

It has recently been shown, in the context of epidemiological models, that dis- 
connection of the interaction network can however be compensated to some 
extent if the striicture of the network itself varies with time, in such a way 
that different parts of the ensemble are not continuously but occasionally in- 
terconnected [4,5,6,7]. Our aim in this paper is to analyze in depth the extreme 
case where, at any given time, each element is connected by the interaction 
network to only one neighbour. Occasional random reconnections can how- 
ever make neighbour to be exchanged. As explained in detail in Section 2 this 
scenario is motivated by the study of a sexually transmitted infection in a 
monogamous population, where each individual has just one sexual partner 
at a time, but partners can change. Specifically, we focus on two critical phe- 
nomena that occur, upon variation of suitable control parameters, in a large, 
connected ensemble of interacting elements, but which are suppressed if the 
interaction network only allows for monogamous static couples. We analyze 
the reappearance of each critical phenomenon when random reconnections are 
allowed to happen. 

In Section 2, we consider the critical transition to endemic states in an epidemi- 
ological model. In Section 3, we study the transition to full synchronization 
in an ensemble of coupled chaotic maps. Both cases are, to a large extent, 
analytically tractable, so that exact results are obtained for the occurrence of 
these critical phenomena in time-varying, highly disconnected networks. 



2 Endemic persistence of SIS epidemics 

The SIS epidemiological model describes spreading of an infection by contagion 
between the members of a population. At any given time, each agent can 
be either susceptible (S) or infectious (I). An S-agent becomes infectious by 
contact with an I-agent, with probability a per time unit. An I-agent, in turn, 
spontaneously becomes susceptible with probability 7 per time unit. Contacts 
between agents are represented by the links of a network. For a static random 
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network -provided that there are no disconnected components, i.e. no portions 
of the population are separated from the rest- the cvohition of the fraction rii 
of I-agents in a large population is well described by the mean field equation 



where ns = 1 — ni is the fraction of S-agents. and ,5 = ka is the infectivity, 
with k the average number of links per agent in the contact network. 

As a function of the infectivity /3, the long-time asymptotic fraction of I- 
agents exhibits two distinct regimes. For ^ < ^, ni vanishes for long times, 
and the infection disappears from the population. For /3 > 7, on the other 
hand, an endemic state with a non-vanishing infection level is reached. The 
asymptotic fraction of infected agents in this state is 



The transition at the critical point where the infectivity /3 equals the recovery 
probability 7 occurs through a transcritical bifurcation [9]. 

Let us now assume that we are dealing with a sexually transmitted infection, 
where contagion can only occur between sexual partners. Suppose also that the 
population is monogamous so that, at any time, each agent has just one part- 
ner. For simplicity, we assume that all agents have partners. In this situation, 
the network of contacts is highly disconnected: for a population comprising 
agents, the network consists of N/2 isolated links defining agent couples. 
If at least one agent is initially infectious within a given couple, the infection 
can persist for some time as a result of repeated contagion between the two 
partners. However, there is a finite probability that in any finite time inter- 
val both agents become susceptible. Assuming that no contacts are allowed 
with other members of the population, the two recovered partners will never 
become infectious again. If couples last forever, consequently, the fate of the 
sexually transmitted infection is to disappear from the population. Due to 
its lack of connectivity, the contact network is unable to sustain a long-time 
endemic state. 

If, on the other hand, sexual partners are allowed to change from time to time, 
even if at any instant the population is still monogamous, the infection may 
spread over the population and, eventually, reach an endemic level. Specifically, 
let us consider that each couple can exchange partners with another 

randomly selected couple {h, I) at a given rate, creating new couples (i, h) and 
(j, /). Would it be possible that, if these recombination events are frequent 
enough, an endemic state is established in the monogamous population? To 
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hi = /3nsni - -fui, 



(1) 




(2) 
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investigate this question we adopt two different approaches, which can to a 
large extent be dealt with analytically. In the first one [7], we compute the 
fraction of I-agents from the contribution of couples formed at different times 
in the past. In the second approach, we study the evolution of the number of 
couples containing two, one, or no infectious partners. 

2.1 Epidemiological dynamics within couples 

Since the moment when two given agents become joined in a couple, the 
probabilities of their being either susceptible or infectious evolve independently 
of the rest of the population. Let po{t), pi{t), and p2{t) be the probabilities 
that, at time t, the couple under study comprises zero, one, and two I-agents, 
respectively. These probabilities satisfy the equations 

Po = IPu 

pi = 27P2 - (7 + (3) 
p2 = - 27P2- 

At all times, po + pi + P2 = 1, so that the analysis can be restricted to the 
system formed by the two last lines of the above equations. As expected, the 
only fixed point of this linear system, (pi,P2) = (0,0), corresponds to the 
disappearance of the infection. The eigenvalues around the fixed point are 

= + 37) ± 1^/32 + 6/37 + 7'. (4) 

Note that < ^+ < 0. The inverse modulus of yU+ gives the typical duration 
time of the infection within a couple, Tj. For 7 ^ /3 and 7 <^ /3 we have, 
respectively, Tj ~ I/7 and Ti ~ P/j'^- The general solution for pi and p2 can 
be written as (pi,P2) = exp(//+i) + A_e_ exp(//_i), with 

e± = f 1, :^^] , (5) 
and where the coefficients A± are determined by the initial values of pi and 

P2- 

When a couple is formed, the initial values of pi and p2 can be evaluated, 
on the basis of mean field-like arguments, from the fraction rii of I-agents 
at that time: pi = 2ni(l — nj) and p2 = nj. Reversing the same arguments, 
the linear combination pi + 2p2 gives the expected fraction of I-agents in the 
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same couple. With these elements, we can now write down an equation for 
the evolution of riiit) taking into account couple recombination events. The 
fraction of I-agents at the present time t is given by the contributions of all 
the present couples, formed at previous times t — s {0 < s < t). The initial 
probabilities within each couple are determined by the fraction of /-agents at 
the time of formation, ni{t — s), and the present probabilities are given by the 
solution to Eqs. (3) discussed above. Summing up all the contributions, we 
get 



ni{t) 



U{t,s) 



2/3 



27 + /X- 



A_[ni{t — s)] exp(/x_s) 



2/3 



27 + )"+, 



A+[ni{t - s)] exp(yLi+s) 



ds. 



(6) 



The coefficients A±{ni) are obtained from the evaluation of the initial proba- 
bilities in a couple in terms of the fraction of I-agents at that time: 



(2ni(l - ni), n^) = A_(ni)e_ + A+{ni)e+. 



(7) 



In Eq. (6), the contribution to ni{t) of the couples formed at t— s is weighted by 
n(t, s), the probability that a couple present at time t has lasted for an interval 
of length s. Since couple recombination occurs at random, this probability 
distribution is a Poissonian function of s, specifically, 

r exp(—rs) , , 

1 — exp(— ri) 



for < s < Here, r is the recombination probability per couple per time 
unit. 

The rather involved form of our equation for ^/(t), Eq. (6), is drastically 
simplified in the long-time limit. We find that for small values of the recom- 
bination rate r, the asymptotic fraction of I-agents vanishes -just as when 
recombination is absent and the infection, confined within couples, eventually 
disappears. There is however a critical value of the recombination rate above 
which an endemic state is reached. The stationary fraction of I-agents is 



The appearance of the endemic state occurs here through a transcritical bi- 
furcation, hke when the infectivity /3 is varied in the standard SIS model. The 
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critical value of the recombination rate is 
27 
^/7-l' 



(10) 



Note that Eq. (9) reduces to Eq. (2) in the limit r — )■ 00. For a given value 
of the infectivity /3, thus, infinitely frequent recombination in a monogamous 
contact pattern sustains a stationary infection level equivalent to that of a 
population with a static, not disconnected pattern. 

To validate the arguments used to derive Eqs. (6) and, in particular, the 
asymptotic result (9), we have performed numerical simulations of the epi- 
demiological model in a recombining monogamous population oi N — 10^ 
agents. Since two couples are involved in each recombination event, the prob- 
ability per time unit of each event in our numerical simulations equals r/2. 
Figure 1 compares analytical and numerical results, which tiun out to show 
excellent agreement. The inset shows the boundary between the regions of 
parameter space where, for asymptotically long times, the infection persists 
or disappears, as given by Eq. (10). 



2.2 Evolution of the number of couples 



Equivalent asymptotic results are obtained from mean field equations for the 
number of couples formed by agents with different epidemiological states. Note 
that, in our problem, the dynamics of the number of couples of different kinds 
corresponds to the evolution of links in the contact network. Let mss, ?7iis, and 
mil be the fraction of couples formed, respectively, by two S-agents, one I- 
agent and one S-agent, and two I-agents. Since at all times mss+^^is+^^ii = 1) 
it is enough to consider the evolution of only two of them, for instance, mss 
and mil. The mean field equations read 

"^ss = \rnl^ - rmssmi + 7"^is, ^^^^ 
rhii = ^mfs - rmssmii + ^mis - 27mii, 



with mis — ^ ~ iTT'SS ~ "^11 ■ The first two terms in the right-hand side of 
both equations stand for the only recombination events that change the num- 
ber of couples of each kind, namely, (SS,11) O (IS, IS). The remaining terms 
correspond to the epidemiological events of contagion and recovery. 

Equations (11) have two stationary solutions. One of them corresponds to the 
disappearance of the infection: 

mss = 1, mil = 0. (12) 
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Fig. 1. Long— time asymptotic fraction of infectious agents function of the 

recombination rate r (normalized by the recovery probability 7), for three values of 
the normalized infectivity (3/'y. Curves and dots correspond, respectively, to analyt- 
ical and numerical results. The inset shows, in the parameter plane (r/7,/3/7), the 
boundary between the regimes of endemic infection and infection disappearance. 

The other fixed point, 

"^ss = ^(2^ + 27 + r), m„ = ^[(^ - l)r - 2^% (13) 

corresponds to the endemic state. In agreement with the results of Section 2.1, 
the two stationary solutions interchange stability when the recombination rate 
attains the critical value given by Eq. (10). The infection level in the endemic 

state can be evaluated from the number of couples as rii = mu + mis/2 = 
(1 + mil — ^ss)/2 and coincides with the fraction of l~agents of Eq. (9). 



3 Synchronization of coupled chaotic maps 

The reappearance of the transition to endemic states in the monogamous SIS 
epidemiological model upon couple recombination events, opens the question 
of whether recombination may compensate disconnection for the occurrence 
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of similar critical phenomena in other ensembles of interacting dynamical ele- 
ments. In this section, we explore this question for the synchronization tran- 
sition of coupled chaotic maps. 

Consider an ensemble of N identical chaotic maps whose individual states 
Xi (i = 1, . . . , N) evolve, in the absence of coupling, according to Xi{t + 1) = 
f[xi{t)]. Let A > 1 be the corresponding Lyapunov coefficient. On the average, 
the distance between two neighboring orbits of a map thus grows by a factor 
A at each iteration. Global (all-to-all) coupling between maps is introduced 
following the standard linear scheme [10,11] 

e ^ 

x,{t + 1) = (1 - e)f[x,{t)] + -Y: f[xj{t)], (14) 

i=i 



where e G [0, 1] is the coupling intensity. It is well known that, if coupling is 
strong enough, the ensemble undergoes a transition to full synchronization. 
Specifically, if 

£ > £e = 1 - ^, (15) 



coupling is able to overcome the exponential separation of chaotic orbits, and 
all the maps converge asymptotically to exactly the same trajectory: Xi{t) = 
Xj{t) for all i, j and t ^ oo [11]. This synchronized trajectory is still chaotic 
-in fact, it reproduces the orbit of a single map- but the motion of the maps 
with respect to each other has been suppressed. Note that the threshold Sc for 
the synchronization transition does not depend on the ensemble size N . 

Consider now a "monogamous" coupling pattern where, at any given time, 
each map is coupled to just one partner. Following the scheme of Eq. (14), if 
maps a and h form a couple, their states evolve according to 

xa{t + 1) = (1 - e)f[xa{t)] + I {f[xa{t)] + f[x,m , ^^^^ 
x,{t + 1) = (1 - e)f[x,{t)] + I {f[xa{t)] + f[x,m . 



If £ > Ec-, the two maps will synchronize, converging asymptotically to the 
same chaotic orbit. Because of the nature of chaotic motion, however, their 
synchronized trajectory will differ from the orbits of other couples. 

Would it be possible that, if couples of maps recombine exchanging partners 

at random at a certain rate, full synchronization all over the ensemble is re- 
covered? By analogy with the occurrence of endemic states in SIS epidemics 
over a recombining monogamous population, we may expect that synchroniza- 
tion spreads over the ensemble if recombination is faster enough, entraining 
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increasingly many maps towards the fully synchronized orbit. However, a dif- 
ference with the case of SIS epidemic is that, now. the recombination rate has 
a limit: one random recombination per couple per iteration step produces the 
maximal possible rearrangement of the coupling pattern. 

3.1 Synchronization dynamics within couples 

To investigate whether full synchronization is possible in the monogamous 
coupling pattern under the effects of recombination, we assume that the en- 
semble is symmetrically concentrated around a reference orbit XQ{t) of the 
map f{x), so that the state of each map differs from xo{t) by a small amount: 
Xi{t) — xo{t) +5i{t). If, as time elapses, the concentration around xo{t) grows 
in such a way that Si{t) — ?> for all i, full synchronization will be achieved. 
Equations (16) imply that the displacements of two coupled maps from the 
reference orbit Xo{t) satisfy the linear equations 



where f'{x) is the derivative of f{x). 

Let us first study the simpler, extreme case where all couples recombine at 

each iteration step. Since new partners are chosen at random, the displace- 
ments Sa{t) and db(t) of the two coupled maps are uncorrelated quantities. In 
other words, the average ({t) = {Sa{t)6h{t)) , performed over all the couples 
in the ensemble, vanishes. Consequently, from Eqs. (16), the variance of the 
displacement over the ensemble, (T^(i) = {^l{t)), is governed by 



Taking into account that the long-time geometric mean value of f'[xo{t)]'^ is 



the variance will asymptotically vanish if (1 — £ -|- 6^/2) < 1 or, equiva- 
lently, if 



Sa{t + 1) = (1 - l)f'[xo{t)]5a{t) + lf[xo{t)]5,{t), 

5,{t + 1) = (1 - l)nxo{t)]6,{t) + lf[xo{t)]5a{t), 



(17) 




(18) 




(19) 



e > El 



= 1 - 



A 



(20) 
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Thus, if the couphng intensity e is larger than the critical value £i, full synchro- 
nization is stable. Note that ei > Ec for any A > 1, so that full synchronization 
for maximal recombination rate requires stronger couplings than when the en- 
semble is globally coupled. Also, there is a critical limit Ai = \/2 ^ 1.41 for the 
Lyapunov coefficient such that, if A > Ai, full synchronization is impossible 
even under the maximal recombination rate. 

When recombination does not occur for all couples at all times, the joint evolu- 
tion of two coupled maps introduces correlations between their displacements 
from the reference orbit. In this case, Eqs. (16) imply that C,{t) — {5a{t)5b{t)) 
and cr'^{t) — {Sl{t)) are governed by 

at +1) = (1 - l)ef'[xo{t)]'a^{t) + {l-e + f )/'K(t)]^CW, ^^1) 
a\t + l) = {l-e+ ^4)r[xo{tWa\t) + (1 - l)ef'[xo{m{t). 



This linear system can be solved exactly. For the couples formed at a certain 
time to, we have C(^o) = and o"^(to) = o'o- With such initial conditions, and 
taking into account Eq. (19), the solution to Eqs. (21) for reads 

a^t) = ^A2(*-*°) [l + (1 - £)2(*-*°)l al (22) 



for i ^ Iq. This quantity is the variance of the displacement from x^ii) of 
those maps whose present couples have formed at time and lasted until the 
present time t. 

Suppose now that, at a certain time r, the ensemble around Xq{t) has variance 
cr^. To evaluate the variance d^^^ at the next time step, we think of the 
ensemble as made up of subensembles consisting of the maps whose present 
couples have formed at times r, r — 1, r — 2, and so on. The fraction of maps in 
the subensemble corresponding to couples formed at time T—n (n = 0, 1, 2, . . .) 
is g„ = r(l — r)", where r is the probability that any couple forms at any given 
time step (0 < r < 1). Using Eq. (22) with t = r and to = ^ — we find that 
the variance of this subensemble, which we denote cr^, changes from o\{t) — 
to 

Assuming that the ensemble has been evolving since an asymptotically long 
time ago, its variance at time r -|- 1 is given by 

-^+1 = E In'^nir + 1) = rX^ E(l - ^r tl' J.,^ -I- (24) 

n=0 n=0 "I" ^) 
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The synchronization threshold is therefore given by the condition 

o° 1 + d _ <r ]2{n+l) 

r^%(^-rr :| =1. (26) 



where Sr is the critical coupling intensity. 

Note that, to obtain this result, wc have used Eq. (22) for t—to = n = 0, 1, 2, . . . 
while, strictly, it holds for t ^ only. Consequently, the threshold condition 
(25) is just an approximation whose validity must be ascertained for each 
specific choice of the map f{x). Below, we present numerical results for a case 
where Eq. (22) holds at any time. 

Whereas, in general, the summation in Eq. (25) cannot be exactly performed, 
two special cases are readily obtained. The first one corresponds to the max- 
imal recombination rate, r = 1, for which we reobtain the critical coupling 
intensity Si of Eq. (20). The second special case gives the recombination rate 
for which the threshold coupling intensity is maximal, Er — 1, namely, 

Tmin = 2 (l - 1) . (26) 



For recombination rates below rmm, synchronization is not possible even for 
the strongest coupling. 



3.2 Numerical simulations 



To test the above results we have performed numerical simulations of recom- 
bining monogamous ensembles of iV = 10'^ chaotic maps for the case of the 
tent map [9]: 



px for < X < 
p{l — x) for I < x < 1, 



(27) 



with X e [0,1] and p e [0,2]. The Lyapunov coefficient of the tent map is 
X — p, which imphes that the dynamics is chaotic for 1 < p < 2. Moreover, 
f'{x)^ — p^ — for all X, so that in Eq. (19) the limit of T ^ oo can be 
dropped, and Eq. (22) holds for any t. 

In our simulations, after joint evolution of the coupled maps following Eq. 
(16), each one of the N/2 couples is allowed to recombine with probabihty p. 
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Fig. 2. Full-synchronization threshold for a recombining monogamous ensemble of 
tent maps, in the parameter plane (p, e), for four values of the map parameter p. Full 
synchronization is stable above and to the right of each boundary. Dots and curves 
stand, respectively, for numerical results on ensembles of 10^ maps and the analytical 
result of Eq. (25). The numerical parameter p is related to the recombination rate 
r through Eq. (28). 

exchanging partners with another randomly selected couple. Since two cou- 
ples are involved at each recombination event, the resulting recombination 
probability per couple per time step is larger than p, and reads 

r = l-(l-p)exp(-p). (28) 



Fully synchronized ensembles are detected by numerically measuring the vari- 
ance crl{t) — N~^J2i[xi{t) — x{t)]'^ with respect to the average state x{t) — 
Si Xi{t). Realizations where cr^(t) falls to the level of numerical round-off 
errors are identified with full synchronization. 

Figure 2 shows the boundary of full synchronization in the parameter space 
(p, e) for four values of the tent map slope p. For each value of p, full synchro- 
nization is stable for coupling intensities e above the boundary. Dots stand 
for the numerical determination of the synchronization threshold and curves 
correspond to the analytical prediction, Eq. (25). The agreement is very good. 
The systematic difference between numerical and analytical results, more vis- 
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Fig. 3. Dispersion of a recombining monogamous ensemble of N = 10^ tent maps 
as a function of the recombination parameter p, for three values ofp and s. Vertical 
dashed lines stand for the numerical determination of the critical value of p, for 
each combination of p and s. Dotted curves are splines added as a guide to the eye. 

ible for small p, is to be ascribed to the numerical discretization of parameter 
space in the search for the synchronization threshold. 

In Fig. 3, we present numerical measurements of the long-time asymptotic 
values of the dispersion function of p, for three combinations of p and 

e. Note that the drop of Ux at the critical recombination rate is abrupt but 
continuous. 



3.3 Two-state symbolic dynamics 



A description of the dynamics of the ensemble of coupled chaotic maps in 
terms of the evolution of different kinds of couple, along the lines developed 
in Section 2.2 for SIS epidemics, is in principle not possible. In fact, while 
the agents in the SIS population admit just two possible individual states - 
either susceptible or infectious, which determine three kinds of couple- the 
individual state of each chaotic map can span a continuous interval. However, 
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as we show in this section, a symbohc representation of each map as a two-state 
variable -either synchronized or unsynchronized- together with a reduced set 
of transition rules between the two states, is able to qualitatively reproduce 
the collective effects of recombination on the synchronization transition. This 
schematic approach has the virtue of capturing the essential mechanisms of 
the interplay between recombination and synchronization, and can also be 
reinterpreted in the context of the SIS epidemics model. 

Let us thus assume that, at any given time, each chaotic map adopts one of 
two states, namely, unsynchronized (U) or synchronized (S). We stress that 
synchronization of an individual map i is here understood as defined with 
respect to the ensemble, for instance, if the variable Xi{t) is within a certain 
small distance from the average x{t) = N^^ J2j Xj{t). In the spirit of Eqs. (11), 
we propose for the evolution of the fractions of couples muu, ^us, and mgs 
the Ansatz 

"^uu = i"^us - rmss^uu + 7"^us - ?7(r)muu, , . 

mss = i"^^s - rmssmuu + /5(r)mus + ?7(r)muu, 



with mus — ^ — ?7T.ss — ''tiuu, and where /5(r), 7, and 77 (r) are non-negative 
quantities. The first two terms in the right-hand side of both equations de- 
scribe the effect of recombination at rate r per couple. They have exactly 
the same origin as in Eqs. (11). The remaining terms, described in detail in 
the following, stand for the transitions induced by the joint evolution of cou- 
pled maps. Note that we are assuming that all these transitions occur always 
toward states where the two maps of a couple are either synchronized or un- 
synchronized. Also, we suppose that the two maps of an SS-couple, as long as 
their link lasts, cannot become unsynchronized from the ensemble. 

In the first place, as time elapses, a US-couple can become either an SS- 
couple or a UU-couple. This describes the tendency towards synchronization 
within each couple. The rate for the transition US — > SS, which we have 
called /3 in Eqs. (29), cannot however be a constant. If it were a constant, in 
fact, all US-couples would become SS-couples for sufficiently long times, thus 
synchronizing with the ensemble even in the absence of recombination. This is 
not the case, though: when recombination is infrequent, maps forming a long- 
lasting couple synchronize to each other, but are generally not synchronized 
to the ensemble. We represent this effect phenomenologically, by ascribing to 
P a dependence on r. The function /3(r) increases as r grows, starting from 
^(0) = 0. 

The transition US — )■ UU, on the other hand, implies the desynchronization 
of a map with respect to the ensemble by interaction with its already unsyn- 
chronized partner. This event does not require recombination: it is rather the 



14 



generally expected outcome within a US couple. We thus assume that its rate, 
7, is constant. 



Finally, Eqs. (29) include, in the last terms of their right-hand side, the tran- 
sition UU SS. Because of the same reasons as in the case of the transition 
US — 7> SS, which also involves the synchronization with the ensemble of previ- 
ously unsynchronized maps, we expect that the corresponding rate, ri{r), is a 
growing function of the recombination rate r. Possibly, however, 77 (r) is much 
smaller that 7(r), because becoming an SS-couple should occur less frequently 
for a UU-couple than for a US-couple. 

Equations (29) have two fixed points. One of them, mss = 1, muu = 0, corre- 
sponds to full synchronization of the ensemble of chaotic maps. For the other 
fixed point, both mss and muu are different from zero, thus corresponding to 
a state of partial synchronization. There is a broad choice of functional forms 
for f3{r) and ri{r) such that these stationary solutions have the expected be- 
haviour -namely, that full synchronization is unstable for small r and becomes 
stable above a critical value of the recombination rate, while the other fixed 
point exhibits opposite stability properties. For the sake of concreteness, let 
us assume the linear dependences /3(r) = (SqT and ri{r) = rj^r. Analysis of the 
eigenvalues of Eqs. (29) around the fixed points shows that, irrespectively of 
the values of /3o and 7, it is enough that rjQ < 1 for the occurrence of a transi- 
tion to full synchronization as r grows. The critical value of the recombination 
rate is 

r.= ^i^. (30) 



Let us now attempt a semi-quantitative comparison of this result with our 
results for the synchronization transition in recombining monogamous ensem- 
bles of chaotic maps, for instance, those depicted in Fig. 2 for tent maps. Since 
7 is the rate for the transition US — )■ UU -where, by interaction within its cou- 
ple, a map becomes desynchronized from the ensemble- it may be identified 
with a measure of the rate at which chaotic orbits diverge from each other, i.e. 
of the Lyapunov coefficient A. The factor ^0, in turn, weights the rate of the 
transition US — >■ SS, where a map becomes synchronized to the ensemble by 
interaction with its already synchronized partner. Thus, /3o is related to the 
strength of the interaction between coupled maps, i.e. to the coupling intensity 
e. The factor r/o plays a similar role but, as mentioned above, its effect should 
be quantitatively less important than that of ^o- Figure 4 is a plot of the 
synchronization threshold as given by Eq. (30). To stress the correspondence 
with Fig. 2, we have plotted (3^ (a measure of the coupling intensity 
function of the recombination rate r (directly related, through Eq. (28), to 
the parameter p) for three values of 7 (a measure of the Lyapunov coefficient, 
\ — piov the tent map), and fixed 770- The semi-quantitative analogy with our 
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Fig. 4. Synchronization thresliold in the parameter plane (r, /3o) for the two-state 
model of coupled chaotic maps, corresponding to three values of the rate 7 and 
r/o = 0.01. Full synchronization is stable in the region of large r and /Sq. Compare 
this plot with Fig. 2. 

numerical results for ensembles of tent maps, shown in Fig. 2, is apparent. 

An analogy between this two-state approach to synchronization and the SIS 
epidemiological model discussed in Section 2 can in turn be derived from the 
identification of the unsynchronized state with the susceptible state on one 
hand, and of the synchronized state and the infectious state on the other. 
With this identification, the transitions US — >■ SS and US — )■ UU correspond, 
respectively, to contagion from the infectious partner and to the spontaneous 
recovery of an infectious agent. The rates /3 and 7 play the same role in both 
dynamical models, with the difference that in the SIS model (3 does not depend 
on the recombination rate. Another difference is that the two-state approach 
to synchronization includes the transition UU — >■ SS, which in the SIS model 
would correspond to simultaneous infection of two susceptible partners -a 
forbidden event. The SIS model, in turn, allows for the recovery of just one 
among two infectious partners, which would stand for the inexistent transition 
SS — >■ US. Notwithstanding these differences, we reahze that the plot of Fig. 4 
is the equivalent in the two-state approach to synchronization as the endemic 
threshold depicted in the insert of Fig. 1 -where both the infectivity /3 and 
the recombination r are normalized by the recovery probability 7. The two 
thresholds have qualitatively the same functional dependence on the respective 
parameter. In addition, the transition to full synchronization in our two-state 
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formulation is a transcritical bifurcation, the same kind of critical phenomenon 
as in the SIS model. 



4 Conclusion 

We have here studied two critical phenomena in the collective behaviour of 
large ensembles of interacting dynamical elements -namely, the appearance of 
endemic states in an SIS epidemiological model, and the stabilization of full 
synchronization of identical chaotic maps- when the corresponding interaction 
patterns are highly disconnected and, concurrently, change with time. Specifi- 
cally, we have considered '"monogamous" interaction patterns, represented by 
networks where each site is connected to only one neighbour at a time, but 
such that neighbours can be exchanged at random at a specified rate. While in 
the absence of neighbour exchange -or, as we have called it, of recombination- 
the occurrence of endemic states and synchronization would be impossible due 
to the lack of connectivity in the ensemble, sufficiently frequent recombination 
events make it possible that coherence is established all over the system, thus 
allowing for organized collective dynamics. 

Recombination of interacting couples in a monogamous pattern introduces a 
new dynamical parameter -the recombination rate. The two critical phenom- 
ena studied here take now place upon variation of this parameter: endemic 
states and full synchronization occur above a certain critical value of the re- 
combination rate. Moreover, for the SIS epidemiological model with a given 
infectivity, we have found that the limit of infinitely large recombination rate 
is equivalent to the situation where the interaction pattern is static but not 
disconnected. For chaotic maps, on the other hand, the fact that time elapses 
by discrete steps imposes a limit to the recombination rate: any interacting 
couple can at most be recombined once per step. This establishes in turn an 
upper limit for the Lyapunov coefficient of individual maps such that they can 
be synchronized by recombination. If the maps are "too chaotic" , synchroniza- 
tion is not possible even at the maximal rate of one recombination per couple 
per time step. 

The fact that, instantaneously, each dynamical element of the two ensembles 
considered here has just only one interaction partner, makes the correspond- 
ing problems analytically tractable to a large extent. In particular, we have 
obtained analytical approximations for the critical values of the recombina- 
tion rate which are in very good agreement with numerical results. For SIS 
epidemics, the critical recombination rate was found both from the integra- 
tion over the whole population of the epidemiological dynamics within each 
couple, and from the dynamics of the number of couples in each epidemiolog- 
ical state. For the chaotic maps, we have replaced this latter approach by a 
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kind of symbolic schematic representation of the dynamics of couples. In this 
representation, which can also be adapted to the SIS model and qualitatively 
reproduces the critical behaviour of the two systems, maps (or epidemiological 
agents) are though of as two-state elements. Their interactions induce transi- 
tions between the two states following a small set of intuitive rules. We con- 
jecture that this kind of representation is capturing the essential mechanisms 
that govern the relative prevalence of unsynchronized/susceptible elements at 
one side of the critical point, and of synchronized/infectious elements at the 
other. 

Collective dynamics on monogamous interaction patterns have the advantage 
of analytical tractability. It should however be borne in mind that these pat- 
terns represent a kind of extreme case among disconnected networks: they have 
the minimal number of links that avoids isolated elements. Networks with less 
severe lack of connectivity should impose lower limitations to the development 
of collective self-organized behaviour. In these cases, therefore, we expect that 
recombination, even at lower rates, will also be able to compensate the lack 
of connectivity, triggering critical phenomena such as those studied here. 
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